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1 Introduction 

The rare decay B — » X s i + £~ offers the B factories the 
possibility to open a new interesting window on flavour 
physics. A measurement of the dilepton spectrum and 
of the forward-backward asymmetry in this decay pro- 
cess has been proposed as a way to test for new physics 
and discriminating between different new physics scenar- 
ios (for a recent review see J]). Indeed, the short-distance- 
dominated flavour-changing neutral-current amplitude of 
B — > X s £ + £~ is extremely sensitive to possible new de- 
grees of freedom, even if these appear well above the elec- 
troweak scale. This type of indirect search for new physics 
relies both on an accurate measurement of the dilepton 
spectrum, and on an accurate theoretical calculation of 
the decay probability. The inclusive B — > X s £ + £~ transi- 
tion just starts to be accessible at the B factories: BELLE 
and BABAR have already data on the rate based on a 
semi-inclusive analysis J3I3E]- These first measurements 
are in agreement with SM expectations, but are still af- 
fected by a 30% error: substantial improvements can be 
expected in the near future. 

From the theoretical point of view, inclusive rare decay 
modes like B — > X s £ + £~ are very attractive because, in 
contrast to most of the exclusive channels, they are theo- 
retically clean observables dominated by the partonic con- 
tributions. Non-perturbative effects in these transitions 
are small and can be systematically accounted for, through 
an expansion in inverse powers of the heavy b quark mass. 
In the specific case of B — > X s £ + £~ , the latter statement is 
applicable only if the cc resonances that show up as large 
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peaks in the dilepton invariant mass spectrum (see fig. ^) 
are removed by appropriate kinematic cuts. In the per- 
turbative windows, namely in the region below and in the 
one above the resonances, theoretical predictions for the 
invariant mass spectrum are dominated by the purely per- 
turbative contributions, and a theoretical precision com- 
parable with the one reached in the decay B — > X s ^ is 
in principle possible 5 ,6 . Regarding the choice of pre- 
cise cuts in the dilepton mass spectrum, it is important to 
stress that the maximal precision is obtained when theory 
and experiments are compared employing the same energy 
cuts and avoiding any kind of extrapolation. 

In these processes QCD corrections lead to a sizable 
modification of the pure short-distance electroweak con- 
tribution, generating large logarithms of the form a"(m(,) 
x log m (m b /Mhoavy), where M hoavy = 0(M W ) and m <n 
(with n = 0, 1, 2, ...). These effects are induced by hard- 
gluon exchange between the quark lines of the one-loop 
electroweak diagrams. The most suitable framework for 
their necessary resummations is an effective low-energy 
theory with five quarks, obtained by integrating out the 
heavy degrees of freedom. Renormalization-group (RG) 
techniques allow for the resummation of the series of lead- 
ing logarithms (LL), a"(ms) log™ (mfe/M), next-to-leading 
logarithms (NLL), a" +1 (m(,) log™ (mj/M), and so on. 



2 NNLL Calculation 

For a detailed discussion of the present status of the per- 
turbative contributions to decay rate and FB asymmetry 
of B — ► X s £ + £~ we refer to []]. Here we simply recall 
that the complete NLL contributions to the decay am- 
plitude can be found in 0. Since the LL contribution to 
the rate turns out to be numerically rather small, NLL 
terms represent an O(l) correction and a computation of 
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Fig. 1. Schematic dilepton mass spectrum of B — > X s £ + £ 
the dashed line corresponds to the perturbative contribution 



NNLL terms is needed if we aim for a numerical accuracy 
below 10%, similar to the one achieved by the NLL calcu- 
lation of B — ► X s ^j. Large parts of the latter can be taken 
over and used in the NNLL calculation of B — > X s £ + £~. 
Thanks to the joint effort of several groups [8*ll§lllUlllll 
I12II13II14II15| . the necessary additional NNLL calculations 
have been practically finalized by now. 

The computation of all the missing initial conditions to 
NNLL precision has been presented in Ref. [BJ. The most 
relevant missing piece of the anomalous dimension ma- 
trix, namely the three-loop mixing of four-quark operators 
into semileptonic ones, has been obtained very recently in 
Ref. 9 . Thanks to these works, the large matching scale 
uncertainty of around 16% of the NLL prediction has been 
removed. 

The two-loop matrix elements of the four-quark oper- 
ators are probably the most difficult part of the NNLL 
enterprise. Asatrian et al. succeeded in using a mass and 
momentum double expansion of the virtual two-loop dia- 
grams |10) . This calculation is based on expansion tech- 
niques which can be applied only in the lower perturbative 
window of the dilepton spectrum. Once the cc threshold 
is reached, the momentum expansion is not a valid proce- 
dure anymore. We have recently extended this calculation 
to the upper perturbative window above the cc thresh- 
old 1141115) : in the next section we shall present the first 
phenomenological outcomes of this analysis. We resorted 
to semi-numerical methods, which are valid both below 
and above the threshold. Regarding the lower perturbative 
window, our analysis provides as an independent confir- 
mation of jlOj , which is particularly welcome in view of its 
complexity: below the cc threshold the agreement between 
our results and those of JU| is excellent. As shown in |10| . 
the NNLL matrix element contributions are a fundamen- 
tal ingredient to reduce the perturbative uncertainty: in 
the lower window the (low) scale dependence gets reduced 
from 13% to 6.5%. As we shall discuss in the next section, 
the residual scale dependence in the upper window is even 
smaller, around the 3% level. 



Another independent ingredient of the NNLL anal- 
ysis is represented by bremsstrahlung corrections (and 
corresponding virtual soft-gluon terms). Also this part 
of the calculation is now available and cross-checked for 
both dilepton spectrum and FB asymmetry |10U11II12II13| . 
While the NNLL program for the FB asymmetry is fully 
completed [12,13 , in principle there are some pieces still 
missing for the integrated dilepton spectrum; however, all 
of them are estimated to be below 1% at the branching 
ratio level. Note that, at this level of accuracy, other sub- 
leading effects become more important. For instance, fur- 
ther studies regarding higher-order electromagnetic effects 
should deem necessary. 



3 Results 

Before discussing the numerical predictions for the inte- 
grated branching ratios, it is worth to emphasize that low- 
and high-dilepton mass regions have complementary vir- 
tues and disadvantages. These can be summarized as fol- 
lows (<j 2 = Mf +e _): 

Virtues of the low-q 2 region: reliable q 2 spectrum; small 
1 /mb corrections; sensitivity to the interference of CV 
and Cg; high rate. 

Disadvantages of the low-q 2 region: difficult to per- 
form a fully inclusive measurement (severe cuts on the 
dilepton energy and/or the hadronic invariant mass); 
long-distance effects due to processes of the type B — ► 
H/X s — > X s + X'£ + £~ not fully under control; non- 
negligible scale and m c dependence. 



Virtues of the high-q 2 region: negligible scale and m c 
dependence due to the strong sensitivity to the Wil- 
son coefficient |Cio| 2 ; easier to perform a fully inclusive 
measurement (small hadronic invariant mass); negli- 
gible long-distance effects of the type B — > fyX s — > 
X s +X'£+£-. 

Disadvantages of the high-q 2 region: q 2 spectrum not 
reliable; sizable l/m& corrections; low rate. 

Given this situation, we believe that future experiments 
should try to measure the branching ratios in both re- 
gions and report separately the two results. These two 
measurements arc indeed affected by different systematic 
uncertainties (of theoretical nature) but they provide dif- 
ferent short-distance information. 

In order to obtain theoretical predictions that can be 
confronted with experiments, it is necessary to convert the 
s = q 2 jm\ range into a range for the measurable dilep- 
ton invariant mass q 2 . Concerning the low-g 2 region, we 
propose as reference interval the range q 2 S [1,6] GeV 2 . 
The lower bound on q 2 is imposed in order to cut a re- 
gion where there is no new information with respect to 
B — > X s ^f and where we cannot trivially combine electron 
and muon modes. Taking into account the input values in 
Tabled the NNLL prediction within the SM is: 
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ml olc = (4.9 ±0.1) GeV 


m c /m b = 0.29 ± 0.02 


a s (M z ) = 0.119 


a em = 1/128 


p = (5.0±|:g) GeV 


1^/^1 = 0.97 



Table 1. Main input values used in the numerical analysis. 



1±. 



L ±6.5% ±2% 

I "si 1 M 1 



±3%! 



rrib (cuts) 



±(4.5±2)%| 1/m , -(1.5±3)%| 
(1.52 ±0.18) x 10~ 5 . 



(1) 



The error denoted by r s \ corresponds to the theoreti- 
cal uncertainty implied by the r(B — > X c ev) normaliza- 
tion which, in turn, is dominated by the uncertainty on 
m c . As already stressed in Ref . [H3 , at present this is the 
dominant source of uncertainty in the low-g 2 region. In 
principle, alternative normalizations such as the one pro- 
posed in Ref. [EJ using B — > X u lv could be used to reduce 
this uncertainty in the future; however, in practice this is 
still the best we can do at the moment. Using the world 
average F(B -> X c ev) = (10.2 ± 0.4)%, we finally obtain: 



BR(q 2 G [1,6] GeV 2 ) = (1.55 ±0.19) x 10" 



(2) 



Concerning the high-dilepton mass region, we propose 
as a reference cut q 2 > 14.4GcV 2 . Using our new NNLL 
evaluation of the two-loop matrix elements (and reanalyz- 
ing nonperturbative effects in this region) we find: 



R, 



high 



dq 



q 2 >14A GeV 2 



4.09 x 10" 



, dr(B X s £+t 
r(B X c ev) 



1±. 



±3% I 



±15%| , t ,-(8±8)%L, ( 2,3)±3%| 

I nib (cuts) v ' I l/m fe I 

(3.76 ±0.72) x 10~ 6 , 



or 



BR{q 2 > 14.4 GeV 2 ) = (3.84 ±0.75) x 10 
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(3) 



(4) 



As can be noted, in this case /j, dependence and intrinsic 
m c dependence induce a negligible uncertainty. In addi- 
tion to the semileptonic normalization (which induces a 
common uncertainty to low- and high-g 2 regions), here 
the largest uncertainties are related to the l/rrif, expan- 
sion. Most of them are parametric, which could be re- 
duced in the future, with the help of more precise data 
on charged-current semileptonic decays B — > X u c lv. The 
leading 15% error denoted by 'mi, (cuts)' indicates the un- 
certainty in the relation between the physical q 2 interval 
and the corresponding interval for the variable s of the 
partonic calculation: this is nothing but the uncertainty 
in the relation between mj and the physical hadron mass. 



(2 3") 

The 8% error denoted by l/rn b ' indicates the combined 
uncertainty due to 1/m 2 and l/w^ corrections: these non- 
perturbative effects induce a divergence in the q 2 spectrum 
(for q 2 — > to 2 ) which cannot be re-absorbed in a shape 
function distribution [6lll7| . This divergence does not pre- 
vent us from making reliable predictions for the integrated 
rate, but it slows down the convergence of the series, which 
turns out to be an effective expansion in inverse powers of 
mb(l — y/s max ), rather than toi,. 1 

We finally note that in both cases (low- and high-g 2 
regions), we have not explicitly indicated the uncertainty 
due to ev cm and |Vi s /Ki,|, which have been fixed to the 
values in Tabled 111 principle, variations of these param- 
eters can be trivially taken into account by appropriate 
multiplicative factors (they both appear as a squared mul- 
tiplicative factors in i? C ut)- However, we stress that a co- 
herent treatment of higher-order electromagnetic effects 
— which is beyond the scope of this work — cannot be 
simply reabsorbed into a redefinition of ct em . 
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For a more detailed discussion of these results we refer the 
reader to |15) . 



